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DIMENSIONAL ANALYSIS USING TORIC IDEALS 

M. A. ATHERTON, R. A. BATES, AND H. P. WYNN 
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» I ' Abstract. Classical dimensional analysis is one of the cornerstones of qual- 

O I itative physics and is also used in the analysis of engineering systems, for 

,^^ ■ example in engineering design. The basic power product relationship in di- 

^ ' mensional analysis is identical to one way of defining toric ideals in algebraic 

^Nj ' geometry, a large and growing field. This paper exploits the toric representa- 

tion to provide a method for automatic dimensional analysis for engineering 
systems. In particular all "primitive", invariants for a particular problem, in 
a well defined sense, can be found using such methods. 

c^ 

-)— > ■ 

C^ ■ 1. Dimensional analysis 

^ \ Dimensional analysis has a long history. It was discussed by Newton and pro- 

O . vided useful intuition to Maxwell, see |9], chapter 3. A recent paper giving a 

jj^ ' pleasant popular overview is [3] . The first rigorous and most well-known treatment 

^~>. is by Buckingham |3], whose name is attached to the main theorem. Dimensional 

1-^ ' analysis is still considered a fundamental part of physics and is taught at an early 

stage in schools and colleges as a basic part of the physics syllabus. It is often 
covered under a heading of qualitative physics [2] • In engineering it gives a useful 
additional tool for the analysis of systems [11] . It is used in engineering design and 
in the formal design of engineering experiments [TU] [H] . It has also been used in 
economics [T]. For an interesting recent application to turbulence and criticality 
\Q\ see [5] [6]. 

VO ' We shall give an algebraic development of dimensional analysis based on the 

T-t^ ' theory of toric ideals and toric varieties. Although this is essentially a reformulation, 

(^ , the algebraic theory itself is by no means elementary. The theory of toric ideals is a 

CO ' live branch of algebraic geometry. We have used [TB] and the recent comprehensive 

volume 8 . We shall see that the methods give all "primitive" invariants for a 
particular problem, in a well-defined sense. 

Within mathematical physics dimensional analysis can also be seen as an ele- 
mentary application of the theory of Lie groups and invariants, when the group is 
C^ , the scale group defined by multiplication. We shall draw on [T5] in the penultimate 

section. 

The basic idea of dimensional analysis is that physical systems use fundamental 
quantities, or units, of mass (M), length (L) and time (T). To this list may been 
added various others such as temperature (K) and current (I), depending on the 
domain. The extent to which new fundamental quantities can be expressed in terms 
of M, L, T goes to the heart of physics but we shall not delve deeply. Mathematical 
models for physical systems use so-called derived quantities such as: force, energy, 
momentum, capacity etc. Dimensional analysis tells us that each one of these 
quantities has units which have a power product representation. Table 1 gives a 
few examples from mechanics. 
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Quantity 


units 


momentum 


MLT-l 


force 


MLT-^ 


work 


ML^T-^ 


energy 


ML^T-^ 


pressure 


ML-^T-^ 


density 


ML-"^ 


volumetric flow 


L^T-^ 



Table 1. Some basic derived quantities 

We note that the formulae for the expression of derived units have integer pow- 
ers. This is critical for our development: it makes them algebraic in the sense of 
polynomial algebra. 

In a physical system we may be interested in a special collection of derived 
quantities. The task of dimensional analysis is to derive dimensionless variables with 
a view to finding, by additional theory or experiment, or by both, the relationship 
between these dimensionless quantities. As mentioned, the key theorem in the area 
is due to Buckingham. In this section we explain it with an example, leaving a 
more detailed discussion until later. 

Rather than use the M, L,T . . . notation we assume that there are some basic 
quantities of interest which we label zi, Z2, ■ ■ ■■ Each quantity is assumed to have 
the scaling property, namely if the fundamental units, which we now call ii, ^2, • • • 
are scaled up or down this induces a transformation on the Zi. Whether this means 
simply a change in units or actual physical scaling of the system is sometimes 
unclear in the literature, but we shall prefer the latter interpretation. 

As example, if zi is force and the fundamental units are mass (ti), length (^2) 
and time (^3), then the scaling transformation is 

With a collection of derived quantities we have one such transformation for each 
Zj. A slightly more realistic formulation is to introduce non-zero constants Cj so 
that in this case we would have 



Zl 



Citit2U Zl, 



but this would make little difference to our derivations 

Here is a well-known example which we shall use as our running example. It 
concerns a body in a fluid and the quantities of interest arc fluid density (zi), 
fluid velocity (22), object diameter (23), fluid viscosity (z4)and fluid resistance (-25). 
Taking the units into account the transformation is: 



/ ^1 \ 

Z2 
(1.1) Z3 

Z4 
\Z5 J 

After a little algebra, or formal use of Buckingham's theorem, we can derive di- 
mensionaless quantities 



t2^tiZ2 
V tit2^t3Z5 J 



2/1 = 2:1 2:2 2:3 2:4 



2/2 



-1 -2 -2 

Z-t ^0 Zo z^ . 
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The first quantity is Reynolds number. The term dimensionless is interpreted by 
saying that replacing each Zj by the j/j in the transformation — s> in 11.11 leaves the 
expression unchanged: the yj are rational invariants of the transformation. The 
dimensionless principal, for our example, embodied in the Buckingham theorem is 
that any function F of xi, . . . , 2:5 which is invariant under -^^ is a function of j/i 
and j/2 which we write: F(yi, ^2)- 

We now sketch the traditional method. The transformation — )■ can be coded up 
by capturing the exponents in the power products. This gives 

/I 1 1 

A= i 0-10-1-2 

\ -3 11-1 1 

This matrix has rank 3 and we can find a full rank 2x5 kernel matrix K. Namely, 
a K which has rank 2 such that A^ K = 0. This is readily computed using existing 
functions in computer algebra such as the "nuUspace" command on Maple. We 
obtained 

1 1 1-10 



K 

' -1 -2 -2 -2 1 

The key point is that the rows of this K give the exponents of 21 , . . . , Z5 in j/i and 
2/2 • However, we can also derive alternative K. For example, 

1 1 1-11 



'0-1 -1 -1 1 

This gives an alternative to y2, above, namely j/s = tj^ t^ t^ t^. The toric approach 
clarifies, among other issues, the immediate problem of the choice of K which this 
example exposes. 

2. TORIC IDEALS 

Algebraic geometry is concerned with ideals and their counterpart algebraic vari- 
eties. We give a very short description here. (Note that we shall use x for variables 
in an abstract algebraic setting reserving z for "real" problems.) A standard ref- 
erence is [8]. We start with the ring of all polynomials in n variables {xi, . . . ,Xn} 
over a field k: k[xi, . . . , a;„]. A set / of polynomials is an ideal if F G / implies 
s{x)f{x) is in / for any s{x) in k[xi, . . . , x„]. By a theorem of Hilbert all ideals are 
finitely generated. That is we can find a set of polynomials fi{x), . . . fm{x) such 
that any f{x) in k[xi, . . . , Xn] can be written f{x) = si{x)fi{x) + • • • + Sm{x)fm{x) 
for some {si{x)} in fc[a;i, . . . ,x„]. An ideal / gives a variety as the set of x such that 
f{x) = for all f{x) S /. The other identity that / is the set of all polynomials 
zero on the variety is not always true, but for the purposes of this paper we will 
alternate freely between varieties and ideal. It will also be enough to work within 
the field Q of rationals. 

Modern computational algebra has benefitted hugely from the theory of Grobner 
bases and the algorithms that grew out of the theory, notably the Buchburger 
algorithm. We will need one more concept, that of a monomial term ordering, or 
term ordering for short. Monomials x" = x"^ ...x"" , where a = ai, . . . , a„ > ie 
tti > 0, i — 1, . . . ,n, drive the theory. A monomial term ordering, written a;" -< x^ 
between is a total (linear) ordering with the addition condition: x" -< x^ implies 
2;a+7 ^ x^'^'^ , for all 7 > 0. Since such an ordering is linear every polynomial / 
has a leading term LT^{f). If we fix the monomial ordering, -<, the Grobner basis 
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G-j = {gi{x), . . . , gmix)} of an ideal / with respect to ^ is a basis such that the 
ideal generated by all monomials in the ideal is the same as that generate by the 
leading terms of G^. Given / and -< the Buchburger algorithm delivers G^. We 
will be concerned with the set of all Grobner bases as -< ranges over all monomial 
term orderings. This is called the fan and is finite, although is can be very large. 

One of the main definitions of a toric ideal fits perfectly with the power product 
transformations of dimensional analysis. It is this observation which motivates this 
paper. We will emphasize the connection by using the same notation: {t,y,A}, 
with X 01 z according to emphasis, but with t, y and A used in both the pure 
algebra and physical theories. 

The following development can be taken from an number of books, but [16^ is 
our main source. The main steps in the definition are. 

(1) The polynomial ring over n variables k[x] = k[xi, . . . , a;„]. 

(2) A d X n matrix A with columns labeled ai, . . . , ad. 

(3) Variables ti, . . . , t^ and the Laurent ring generated by the U and the inverses 
t^^. We write this as 

k[t,t-^]=k[ti,...,td,t^\...t-^]. 

(4) A power product mapping from k[x] to fc[t,t~^] defined by A: 

x^ -> f^' 

The kernel of the mapping in (4) above is the toric ideal. It can be considered 
as the ideal obtained by formally eliminating t t from the ideal: 

(xj - t*^',? = 1, . . .,n) 

The following paragraph should be considered as a theorem. 
The generators of the toric ideal ideal are related to the kernel of A in the follow 
way. The generators are all so-called binomials 

x" - x^, 

where u and v are non-negative integer vectors with the property that 

Au = A-v. 

The last equation can be written A(u — v) = 0, which is equivalent to u — v 
being in the kernel of A. 

The connection with dimensional analysis should now be clear. Let us put di- 
mensional analysis on a similar notational footing, only using z instead of x. Start 
with a d X n matrix A with columns {a^}. The general form of the mapping — )■ in 
II. II becomes 



(2.1) Zi ^t'^'z,, i = l,...,n 
We can write this in matrix terms as 

(2.2) z -^ t^z 

Now, suppose we have a possible invariant yj. Using u, v to denote integer vectors 
with non-negative entries to distinguish the positive from the negative exponents 
and write 

2/j =z"Jz-^J 
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The condition to be an invariant is that substituting each Zj by tjj in the right 
hand side of 12.21 for z leaves yj unchanged. But the condition for this is 

z"z- = (f4z)"(f4z)-\ J = l,...,d, 

which is equivalent to 

Auj ~ A-Vj =0, J = l,...,d, 

exactly the toric condition. We have proved our main result: 

Theorem 2.1. A variable y is a dimensional invariant in a system defined by a 
matrix A, with derived variable z, if any only if it takes the form 

y = z"z-- 

where u and v are non negative integer vectors such that Au — Av. Moreover the 
set of all quantities 

Z — Z , 

is the toric ideal I a with generator matrix A. 

A brief summary is to say that the set of all dimensional quantities y associated 
with A are exactly those given by the toric ideal Ia- 

We can give a minimal set of generators for the toric ideal of our running example. 
We use the "Toric" function on the computer algebra package CoCoa [7J, which 
takes the matrix A as input. Simply to ease the notation in the use of computer 
algebra we use a, . . . , e, for zi, . . . , Z5. The script with output is. 

Use R ::= QQ[a, b, c, d, e]; 
Toric([[l,0,0,l,l], [-3,1,1,-1,1], [0,-1,0,-1,-2]]); 
Ideal(— d^ -I- ae, abc — d, bed — e) 

By the theorem, given any generator we have a invariant. Thus —d^ + ae yields 
^. Thus we have gives three invariants: 

ae abc bed 

We see that the second two ideal generators give exactly the dimensional variables 
from the kernel matrix K', above. A key point is that the toric ideal may have 
more generators than the rank of the kernel in Buckingham's theorem. The next 
section explains why this is so. 

2.1. Saturation and Grobner bases. To summarise, the toric version of dimen- 
sional analysis says that we can generate dimensionless quantities from the toric 
ideal which is the elimination ideal of the original power product representation, 
being careful to use elimination in the proper algebraic sense. 

A lattice ideal associated with an integer defining matrix A is the ideal based on 
a full rank kernel matrix. That is if A is d x n with rank d then we find an integer 
n X d — n matrix K, with rank n ~ d with rows ki, . . . , k„_rf with A'^ K = 0. 

The corresponding lattice ideal is generated by {t'^^ }. For our first if in[T] lattice 
ideal has two generators: 

/ -1 -1 -2 -2 -2 \ 

{ZlZ2Z3Z^ ,Z^ Z^ Zg 2:4 25)- 

But, as we have seen, this has one fewer generators than the toric ideal. However, 
given any such lattice ideal we can obtain the toric ideal using a process called 
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saturation. The process has two steps. Fix the defining matrix A and let Ia be a 
lattice ideal associated with A. 

(1) Select a dummy variable s and adjoin to the lattice ideal the generator 
s YVi=i Xj + 1. That is form the union 

n 

(2) Eliminate s from I'^ to give the toric ideal for {xi, . . . ,a;„}. That is, the 
toric ideal is obtained as the elimination ideal for {a;i, . . . , a;„}. 

The process of elimination in this saturation process is a formal procedure and leads 
to a reduced Grobner basis of the toric ideal which in general depends in general on 
the monomial ordering used in the elimination algorithm. 

This process gives an explanation for the fact that the toric ideal contains, but 
is not necessarily equal to the lattice ideal. Recall that unions of ideals is mirrored 
by intersections of varieties. The addition condition s J^ ■ j^ Xj + 1 = giving the 
variety defined by I*^ forces all the Xj to be nonzero. This property is inherited by 
the toric ideal. It implies that if any Xj — then all Xj and zero. That is to say, 
saturation removes the principal axes and all axial subspaces. 

This gives a nice physical interpretation. If we exclude the origin, then for the 
toric variety associated with the toric ideal must not contain any other zeros. Trans- 
lated into the original Zj variables, the toric ideal description of the dimensionless 
quantities is appropriate when non of the defining variables Zj is allowed to be zero. 
This removal of zeros is intimately connected with the abstract definitions of toric 
varieties based on the concept of a torus in complex variables, but we do not develop 
this here, see [S]. 

2.2. The Grobner fan, primitive invariants and the Graver basis. A natural 
question given the ease of computing invariants using toric methods is whether the 
invariants obtained in this way are in some sense minimal. This turns out to be 
the case. We can illustrate this with our example. A little inspection of the basis 
— d^ + ae, abc — d, bed — e shows that we cannot get simpler invariants from this 
basis by multiplication (or division) : eg if 

ae bed 

then 2/12/2 — ^ which, although a new invariant, is not obtained by reducing the 
numerator or denominator of any of the original invariants. 

Definition 2.2. A basis element z" — z^ of I a is called is called primitive if there 
is no basis element invariant z" — z"" such that such that z" divides z" and z"" 
divides z". We call an invariant y = z^z^'" primitive if and only if z" — z^ is 
primitive as a basis element of I a • 

Lemma 4.6 of [TB] is 

Theorem 2.3. Every invariant obtained from as a reduced Grobner basis of Ia is 
primitive. 

Note that in what follows we are a little lazy in not disassociating an invariant 
from its inverse. 
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As mentioned, as we range over all monomial term orderings defining the in- 
dividual Grobner basis we obtain the complete Grobner fan and by the lemma 
and our definition all resulting primitive invariants are primitive. This union of 
bases is called the universal Grobner basis and the computer programme Gfan is 
recommended to compute the fan |13) . 

We return to our running example. If we put the G-basis element (— d^+ae, abc— 
d, bed — e) into Gfan we obtain the full fan as 

{bed — e, ae — d^, abc — d), (e — bed, abe ~ d), {d^ — ae, bed — d, abe — d) 

{d — abe, ab c — e), {e — ab e ,d — abe), 

the first of which is the input basis. 

The universal Grobner basis of distinct basis terms (ignoring the sign change) 
copy is: 

bed — e,ae — d , abc — d,ab c — e 
and we have a new primitive invariant: '^ . 



The set of all primitive polynomials, which may be larger than that giving the 
imion of the basis elements in the fan, is called the Graver basis. Algorithm 7.2 of 
[l6 | can be used for this. 

Briefly, the method starts by constructing from A an extended matrix called the 
Lawrence lifting: 

' A Q ' 
I I ' 



A 



where the zero is a d x n zero matrix and / is a d x d identity matrix. Then 
introducing n more derived variables to make a set zi, . . . , Zn, Zn+i, • ■ • 22n a toric 
ideal is constructed using A. Finally, set Zn+i = ■ ■ ■ = Z2n — 1 

The method is conveniently set out in the help screen of "ToricIdealBasis" on 
Maple. After inputting 



A 



1 








1 


1 




















-1 





-1 


-2 

















3 


1 


1 


-1 


1 

















1 














1 

















1 














1 

















1 














1 

















1 














1 

















1 














1 



we use the commands 



zs :— [seq{z[i],i — 1..10)]; 

T := TorieIdealBasis(A, zs, plex{op{zs)) , method =' hs' , grading = grd); 

G := subs{[seq{zs[i\ = l,i = 6..10)],T); 



This yields 



(Z2Z3Z4 - Z5,ziZ5 - Z4,-2;4 + ziZ22;3,-Z5 + 2122^3) 



In this case the set is the same as given by the fan. That is, the universsal Grobner 
basis and the Graver basis are the same. 
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3. Further examples 

For each of the examples we give the derived quantities using the classical nota- 
tion, (i) the A matrix (ii) a single toric ideal basis give by the default function on 
CoCoa and (iii) a full set of primitive basis elements, that is the Graver basis, given 
by the maple 'ToricIdealBasis" command. From this a full set of primitive invari- 
ants is immediate. It turns out that for all except one of our examples (windmill) 
the Graver basis is also the universal Grobner basis. We try to mention when we 
find well-known invariants. 



3.1. Windmill. This standard problem is taken from [TT", Section 9.3. (We have 
changed d there to D). 

It concerns a simple windmill widely used to pump water. The table is 



shaft power, P 


ML'Y-'^ 


diameter, D 


L 


wind speed, V 


LT-^ 


rotational speed, n 


T-i 


air density p 


AfL-3 



1 








1 


2 1 


1 





-3 


-3 


-1 


-1 






The A-matrix is 



In the a,b . . . notation we obtain, from Cocoa, a basis with 4 terms : 

{bd — c, 6 c e — a, c e — ad ,bc e — ad) 
The first entry give a dimensionless quantities discussed in the book: 

V 

The universal Grobner basis obtained from the Gfan gives five terms 

{bd — c,b c e ~ a,c e — ad ,bc e — ad, b^d e — a) 

The last of these is also discussed in the book; it gives the invariant 

P 



pn^d^ 
A full set of 7 primitive invariants, the Graver basis, is 

{bd — c, b^c^e — a, C'e — ad^ , bc^e — ad, b^d^e — a, b'^cd^e — a, be e — ad). 

Since rank(A) = 3 there are only two algebraically independent invariants. The 
standard argument may suggest testing the relationship between any two indepen- 
dent invariants, for example in a wind tunnel. An important question, which should 
be the subject of further research, is say which two or, more generally, whether the 
dimensional analysis is sufficiently trusted to test only one pair and infer other 
relationships from the algebra. 
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3.2. Forced convection. The interest is in the following derived quantities: the 
forced convection coefficient h, the velocity, u, the characteristic length of the heat 
transfer surface L, the conductivity of the fluid fc, the viscosity, /i, the fluid specific 
heat capacity, c and the fluid density, p. The fundamental dimensions are M, L, T 
and two new ones temperature [K) and energy (J). With columns in the order of 
the listed the rows in the units order the A-matrix is 



/ 



V 

Resorting to the a, 6, . . . notation we have from, Cocoa, the ideal 















1 


-1 


1 \ 


-2 


1 


1 


-1 


-1 





-3 


-1 


-1 





-1 


-1 








-1 








-1 





-1 





1 








1 





1 


0/ 



{ac — d,ef — d, beg — e, bfg — a,ae — bdg), 

giving invariants: 

ac ef beg bfg ae 
d ^ d ' e ' a ' bdg 

The first three of these are well-known invariants: 



(3.1) 

(3.2) 
(3.3) 



Reynolds number : R 

Nussclt number, N 
Prandtl number, P 



puL 

M 
hL 

~k 
pc 



In preparing this paper it was pleasing to obtain these directly from the computer 
on the first run. The full set of 7 primitive basis elements is 

(ac — d, ef — d, beg — e, bfg — a,ae — bdg, befg — d,ac — ef) 
The simplest of the "new" primitive invariants is from ae — ef: 

hL 
pc ' 

which is the Reynolds/Nusselt. 

3.3. Electrodynamics. As an exercise we take six basic quantities for electro- 
dynamics, used the literature to give some expression in terms of mass (M), length 
(i). Time (L) and current (A). We do not have any particular elctromagnetic 
device in mind, but simply try to find some dimensionless quantities. The table 
below gives one version: 



Quantity 


units 


charge 


TA 


potential 


ML^T-^A-^ 


capacitance 


M-^L-^T^A^ 


inductance 


ML^T-^A-^ 


resistance 


ML^T-^A-^ 
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The A-matrix is 



A 



( 


1 


-1 


1 


1 \ 





2 


-2 


2 


2 


1 


-3 


4 


-2 


-3 


V 1 


-1 


2 


-2 


-2/ 



CoCoa gives 

{be — a, — ce^ + d, — ae^ + 6d). 
Note that A only has rank 3. It turns out that this is a complete list of primitive 
basis elements. 

3.4. Quantum. Toric ideals are embedded in advanced models in physics but one 
can get some way with simple dimensional analysis. This example is given in some 
form by a number of authors. We found [T3], section 1.3.1, useful. The hydrogen 
atom consists of a proton and a neutron and the Bohr radius is the distance between 
them. We have used slightly non-standard notation. In a somewhat cavalier manner 
we have introduced the speed of light as derived quantity. 



mass of electron, me 




M 


Bohr radius, ao 




L 


energy, E 




ML'^T-^ 


Plank's constant, h 




ML^T-^ 


permitivity of vacuum 


(squared), e^ 


ML^T-^ 


speed of light, c 




LT-^ 




Then the ^-matrix is 

1110 
2 2 3 1 
-2 -1 -2 -1 
The ideal is 

{—df + e,bc — df, ate — d , —cd + ae , abf — d,af — c, aef — cd), 

(the algebraic e is e^ and the algebraic c should not be confused with the speed of 
light). The first terms gives an invariant called the "fine structure constant" 

he 

2 . . . 

If we take the third term and interpret ^— being invariant as stating that v = 

2 

constant x ^—, then we have a well known formula for gq interpreted as the size of 
the hydrogen atom: 

flQ = constant x ^. 

TOeC 

We cannot resist stating that the sixth basis element, a/^ — c gives 

E = constant x rrieC . 
The Graver basis gives a full set of 10 primitive invariants for the hydrogen atom 



{—df+e,bc—df,abe—d ,—cd +ae ,abf—d,af —c,aef—cd,bc—e,abf — e, — / +ab c) 
It is not known whether this list has been given explicitly before. 
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4. Group invariance 



Dimensional analysis should be considered as a special case of the theory of 
groups invariance and in an attempt to suggest a natural generalisation we very 
briefly sketch the theory of invariants. 

We start with the action of a Lie group G acting on a manifold M in R'^. 
The manifold will be our model and the group something to do with our physical 
understanding of the physics being modelled. The orbit of 0{x) is a point a; in M 
be the set of all g{x) for all g in G. If M is invariant under G then 0{x) C M . 
This sets up an equivalence relations with members of A/ in the same orbit being 
equivalent. The collection of equivalence classes is denoted by the quotient M/G 
and the projection tt : M — > M/G maps every member of of M into its correct 
equivalence class. If we are lucky then M/G is a manifold in its own right and 
we say that G acts regularly on M . Also, the mapping tt can be used to set up a 
coordinate system on M/G and note that tt itself is an invariant. This discussion 
leads naturally to the following 

Proposition 4.1. Let a group G act regularly on a manifold M. Consider a 
manifold defined by a smooth function F is a set Sp = {x\f{x) = 0. Its is G- 
invariant is and only there is a function F* defining a smooth sub-manifold Sp* = 
{y\F*{y) = 0} on M/G such that 

Sp, — t:{Sp), 

where tt is the projection form hi to M /G . 

A one parameter Lie group G shifts a point x along an integral curve ^(e, a;) 
called a flow If we expand ^(e, x) in a Taylor expansion in e we obtain: 

*(e) = a; + ee(a;)+0(e2). 

The term ^(x) = (^i(a;), . . . , ^d(x) defines a vector field and we can write v in local 
coordinates in classical 

ox\ oxd 

A function ■0 is an invariant if vijj = or 

^'<^'£ --«*)£-■ 

This is a first order partial differential equation which can be solved by writing 

down 

dxi dxi 

£.i{x) S.d{x)' 

namely by the methods of characteristics. The solutions take the form: 

i>l{x) = Ci,--- ^Ipmix) = Cm, 

where the 4'j ^^e the invariants. 

In our notation e becomes t and the mapping — ^ in (jl.l|) is 

*(t,z) =t^z. 

Matrix partial differentiation with respect to t, and setting all ti = 1 gives the 
infinitesmal generators: 

v = A— 
dz 
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An interpretation of the toric variety is as characterising the orbits of the group, 
as discussed above. We have not formally proved the Buckingham theorem, but 
drawing on the above discussion it is give as Theorem 2.22 in |14) . 

5. Discussion 

We have seen that the toric ideal method, via the Graver basis, is a fast way to 
compute all primitive invariants in dimensional analysis. There are three areas of 
further study which this suggests. 

The first area arises from the possibility that different physical systems may yield 
different types of toric ideal or variety. The most important general class is normal 
toric varieties. Briefly such varieties are related to polyhedral cones and polyhedra 
with integer or rational generators. The standard approach is to take the such a 
cone a and compute its Hilbert basis, which is a set of integer generators of the 
dual cone which gives all integer grid points in that cone. From this there is a 
natural toric ideal. But an open problem, it seems to the authors, is whether the 
rich theory of normal and polyhedra has a role in classical physics and engineering. 

The second area would be the natural development from the last section. A 
discussion missing from in this paper is the way in which differentials are convert 
to derived quantities. For example velocity, which is -S, for some length variable 
y and time t is awarded the derive quantity LT^^. One way to keep the advan- 
tages of awarding derived quantities to differential terms, but keep the meaning of 
differentials is to use combinations of differential and polynomial operators. The 
algebraic environment which combines differential operators of this kind with poly- 
nomial algebras are differential algebras and in particular Weyl and Orr algebras. 
It would be useful to develop a type of generalization of dimensional analysis which 
combined differential algebras with the invariance touched on in the last section. 

The third area is considered because the authors came to this work from the use 
of experimental design methods in engineering. It seems that having easy access 
to all primitive invariants should expand the scope of experimental design methods 
based on invariants, which is a small but established field see [11] [12]. This is 
mention briefly in ()3.1|) . The authors hope to develop this idea. 

References 

1. Barnett, W., Dimensions and economics: some problems, Quartertly J. Austrian Economics 
1 (2007), no. 95-104. 

2. Bhaskar, R. and Nigan, A., Qualitative physics using dimensional analysis, Artifivcal Intelli- 
gence 45 (1990), 73-111. 

3. Hershberger, R.E. Bolster, D. and Donelly, R. J., Dynamic silimarity, the dimensionaless 
science, Physices today (2011), 42-46. 

4. Buckingham, E., On physically similar systems: illustration of the use of dimensional anal- 
ysis, Phys, Rev 4 (1914), 345-376. 

5. Chapman, S. C. and Watkins, N. C, Avalanching systems under intermediate driving rate. 
Plasma Phys. Control Fusion (2009), 1-9. 

6. Rowlands, C Chapman, S. C. and JOURNAL = Physics of Plasmas YEAR = 2009 volume = 
16 pages = 012303 Watkins, N. C, TITLE = Macroscopic for control parameter for avalanche 
models for bursty transport . 

7. CoCoATeam, CoCoA; a system for doi ng Computations in Commutative Algebra, Available 
at http://cocoa.dima.uiiige.it. 

8. Little, J. B. Cox, D. A. and Schenck, H.K., Toric varieties, 124 ed.. Graduate Studies in 
Mathematics, American Mathematicval Society, Providence, Rhode Island, 2011. 

9. D'Agostino, S., A history of the ideas of theoretical physics, Kluwer, Dordrecht, 2002. 



DIMENSIONAL ANALYSIS 13 

10. Gibbings, J. C, The systematic experiment, Cambridge University Press, Cambridge, 1986. 

11. , Dimensional analysis, Springer, London, 2011. 

12. Grove, D. M. and Davis, T. P., Engineering, quality and experimental design, Longmans, 
London, 1992. 

13. Jensen, A. N., Gfan, a software system for groebner fans, http;//wwvir.math.tu- 
bcrlin.de/ jensen/software/gfan/gfan.html (2012). 

14. Olver, P. J., Application of lie groups to differential equations. Springer- Verlag, New York, 
1986. 

15. Smith, H., An introduction to quantum mechanics. World Scientific, New York, 1991. 

16. Sturmfels, B., Grobner bases and convex polytopes, 8 ed.. University Lecture Series, American 
Mathematicval Society, Providence, Rhode Island, 1996. 

BauNEL University, Middlesex, UK 
E-mail address: Mark.Athertonabrunel.ac.uk 

Rolls-Royce plc, Derby, UKJ 

E-mail address: Ron.BatesORolls-Royce.com 

London School of Economics, London, UK 
E-mail address: h. wynn01se.ac.uk 



